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Technical Abstract 

We numerically simulate the effects of noise-induced sampling of alternative Hamiltonian 
paths on the ability of quantum adiabatic search (QuAdS) to solve randomly-generated 
instances of the NP-Complete problem N-bit Exact Cover 3. The noise-averaged median 
runtime is determined as noise-power and number of bits N are varied, and power-law and 
exponential fits are made to the data. Noise is seen to slowdown QuAdS, though a downward 
shift in the scaling exponent is found for N > 12 over a range of noise-power values. We 
discuss whether this shift might be connected to arguments in the literature that suggest 
that altering the Hamiltonian path might benefit QuAdS performance. 

Non- Technical Abstract 

We numerically simulate the effects of noise on the ability of the Quantum Adiabatic Search 
(QuAdS) algorithm to solve randomly-generated instances of the NP-Complete problem N- 
bit Exact Cover 3. The noise-averaged median runtime is determined as noise-power and 
number of bits N are varied, and power-law and exponential fits are made to the data. 
Noise is seen to slowdown QuAdS, though a downward shift in the scaling exponent is found 
for N > 12 over a range of noise-power values. We discuss whether this shift might be 
connected to arguments in the literature that suggest that altering the path followed by the 
search Hamiltonian might benefit QuAdS performance. 
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NP-Completeness; noise 



1 Introduction 



One of the deepest open questions in theoretical computer science is whether the computa- 
tional complexity classes P and NP are equal pp. It is widely conjectured that these two 
classes are different, though it is known that should a polynomial-time algorithm be found for 
an ./VP-Complete problem, then P = NP. In 2001, Farhi et. al. [2] examined whether a quan- 
tum algorithm might be able to solve an AP-Complete problem in polynomial-time. They 
used the quantum adiabatic search (QuAdS) algorithm [3] to find solutions to randomly- 
generated hard instances of the NP-Complete problem iV-bit Exact Cover 3 which they 
believed to be classically intractable for sufficiently large N. Using a digital computer they 
numerically simulated the QuAdS dynamics. They determined the algorithm's median run- 
time to solve this class of instances for a restricted range of N values and found their results 
could be fit with a quadratic scaling relation. It was noted that should classical algorithms 
truly require exponential time to solve this class of instances, and should the quadratic scal- 
ing behavior of QuAdS persist to large N, then QuAdS could outperform classical algorithms 
on this class of instances, though not necessarily on the worst case instances. 

QuAdS works well so long as the quantum dynamics is adiabatic. This requires the 
runtime T to be large compared to 1/ A 2 , where A is the smallest value (encountered during 
the dynamical evolution) of the energy gap between the ground and first-excited states. 
When A is too small, the adiabatic condition is violated, and QuAdS performance suffers. 
Farhi and co-workers [I] examined a case where a failure of QuAdS was transformed into a 
success if the path followed by the search Hamiltonian H(t) differed from the Hamiltonian 
path used in Refs. [2] and [3] that linearly interpolates from an initial to a final Hamiltonian. 
The essential point is that, should the linearly interpolating search Hamiltonian produce a 
A that is too small, varying the Hamiltonian path may cause the new search Hamiltonian 
to produce a larger A and thus improve QuAdS performance. 

A number of papers have considered the robustness of QuAdS performance to noise [5] 
[8]. In this paper we extend the simulations reported in Ref. [8] in two important ways. 
First, the simulation results presented here examine QuAdS performance in the presence of 
non-uniform noise in which each qubit interacts with a different noise field. Ref. [5] focused 
on uniform noise. Second, the simulations in this paper are done at larger noise power and 
larger numbers of qubits N. These differences allow the simulations in this paper to sample 
a larger range of Hamiltonian path variations than was possible in Ref. [8] , and so provide a 
better opportunity to explore how Hamiltonian path variation impacts QuAdS performance. 

This paper is organized as follows: (i) Section [2] briefly summarizes the QuAdS algorithm, 
our noise model, and the simulation protocol; (ii) Section [3] presents our simulation results; 
and (iii) Section 0] closes with a discussion of these results. 

2 Background 

We begin with a description of the AP-Complete problem A-bit Exact Cover 3 (EC3) [H[2]. 
An instance of EC3 is specified by a set of clauses C*, with % = 1, . . . , L, and each clause 
Ci is specified by 3 integers: Cj = (a(i),b(i),c(i)). The integers a(i), b(i), and c(i) satisfy 
a{i) < b(i) < c(i) and take values in the range [1, . . . , N]. Generally, the number of clauses 
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L varies from one EC3 instance to another. A binary vector z = (z±, . . . , zn) (with Zj = 0, 1) 
satisfies the clause if its components z a ^, z b ^, and z c (i) satisfy z a ^ + z^i) + £ c (i) — 1- 
Otherwise, z is said to violate Cj. A binary vector z solves an instance of EC3 if it satisfies 
all of its clauses. Finally, an EC3 instance is said to have a unique satisfying assignment 
(USA instance) if only one binary vector solves it. 

In QuAdS an n-qubit register is initially prepared in the groundstate of a Hamiltonian 
Hi. The only conditions placed on Hi are that its groundstate be non-degenerate and easy 
to prepare. Hi is then adiabatically evolved over a time T into a final Hamiltonian Hp. The 
final Hamiltonian is constructed so that a basis for its groundstate eigenspace encodes all 
solutions to the instance of the computational problem that is to be solved. The details of 
how Hp is constructed from an instance of EC3 are described in Ref . [2] . The construction 
of Hi and Hp are such that both are dimensionless and have energy- level spacing AE ~ 1. 
Since the initial state is the groundstate of Hi, the adiabatic dynamics insures that the final 
state will be in the groundstate eigenspace of Hp with probability P — > 1 as T — > oo. An 
appropriate measurement of the quantum register at the end of the adiabatic evolution then 
yields one of the instance solutions. In Refs. [2] and [3] the time-dependent Hamiltonian 
Hit) that drives the QuAdS dynamics linearly interpolates from Hi to Hp, 

H(t) ={ 1 ~f) H i + (f) H p » 

where < t < T, and T is sufficiently large that H(t) produces adiabatic dynamics. The 
simulations in Ref. [2] randomly generated 75 USA instances of iV-bit EC3 which were 
believed to be hard instances for both classical algorithms and QuAdS. The median runtime 
T(N) for QuAdS to succeed on these instances was found for 7 < N < 20. It was found 
that the simulation results could be fit with a quadratic scaling relation T(N) ~ N 2 . 

To study the impact of noise on QuAdS we introduce classical noise fields Nj(t) (J = 
1, . . . , N) that couple to the qubits via the Zeeman interaction 

N 

H int (t) = -J2<Tj-' N A t ) ■ (2) 

In this paper we focus on non-uniform noise where each qubit is acted on by a different 
noise field: Nj(t) ^ Nj(t) (j ^ i). A detailed presentation of our noise model is given in 
Ref. [H] — we summarize its essential properties here. Each noise field Nj(t) is a sequence of 
randomly occurring fluctuations with profile Fj(t): 

N f 

N i (t) = ^F i (t-t fc ) (j = l,...,iV) . (3) 

k=i 

Here tf. is the temporal center of the /c-th fluctuation and Nf is the number of fluctuations. 
The fluctuations have the following statistical properties: (1) the number of fluctuations 
Nf that occur in a time T is Poisson distributed with average fluctuation rate n; (2) each 
fluctuation profile Fj(t — tk) is a square pulse with height Xj^ and temporal width 2r, where 
t is the thermal relaxation time; (3) the height Xj^ is Gaussian distributed with zero mean 
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and variance a 2 ; and (4) the times are uniformly distributed over [0, T\. The simulations 
allow the polarization of Fj(t) to be either: (i) fixed along x, y, or z; or (ii) to fluctuate 
simultaneously along all 3 directions. In Ref. [8] it was found that noise polarized along 
y caused the largest slowdown of QuAdS and so we focus on y-polarized noise throughout 
this paper. The time-averaged noise power P is related [8j to n, a 2 , and r via P = 2na 2 r. 
The simulations described below use a = 0.2; r = 1; and average noise power in the range 
0.001 < P < 0.013. The average fluctuation rate is then determined from n = P / 2a 2 r. 

The QuAdS simulation protocol with noise is described in Ref. [8]. As with noiseless 
QuAdS, the protocol with noise begins by producing 75 randomly generated USA instances 
of A-bit EC3. The simulations described below were done for 7 < N < 16. For each USA 
instance, 10 noise environments {N™"(t) : j = 1, . . . , N; m = 1, . . . , 10} were generated. For 
each m, H int (t) is determined by plugging the {Nj m (t)} into eq. (J2j). The noiseless QuAdS 
Hamiltonian H(t) is given by eq. ([I]), where Hi is the same for all USA instances and each 
USA instance determines its own Hp [2J. The noisy QuAdS Hamiltonian T~C(t) is then 

H(t) = H(t) + H mt (t) ■ (4) 

For each USA instance and noise environment, 7~C{t) drives the Schrodinger dynamics of 
QuAdS. This dynamics is numerically simulated to find the runtime for QuAdS to succeed 
on that instance and that noise environment. For each N, a total of 750 = 75 x 10 QuAdS 
runtimes are generated. The noise-averaged median runtime (T(N)) is then identified with 
the median of the 750 runtimes. We determined the best power-law and exponential fits to 
the simulation results, and calculate their associated xju anc ^ probability P(x 2 > xju)- The 
simulations were done on the National Science Foundation TeraGrid cluster. 

3 Results 

We now present our simulation results. As a baseline for the noisy simulations we repeat 
the noiseless calculation of Ref. [2J for 7 < N < 19. Our results appear in Figure [T] which 
contains best power-law and exponential fits to the data. The power-law fit T(N) = aN b 
has fit parameters a = 0.1016 and b = 2.079. The value of \ 2 f° r the fit is xju = 0.321 
and the probability P(x 2 > xju) = 0.9999. The closer this probability is to 1, the more 
consistent the data-set is with the fitting function. The exponential fit T(N) = a[exp(6A) — 1] 
has: (i) fit parameters a = 4.707 and b = 0.1282; (ii) xju = 0.296; and (iii) probability 
P{x 2 > xju) = 0.9999. Both fits are excellent and the power-law fit is consistent with the 
result of Ref. [2]. 

The results for our noisy QuAdS simulations appear in 3 figures and 2 tables. Figures [2} 
|3] and H] plot the noise-averaged median runtime (T(N)) versus N for average noise power 
P= (0.001,0.003); (0.00_5, 0.007); and (0.009,0.013), respectively. As in Figured each plot 
contains a power-law ((T(N)) = aN b ) and exponential ((T(N)) = a[exp(6A) — 1]) fit to 
the simulation results. The parameters associated with the power-law and exponential fits 
appear, respectively, in Tables [T] and [2J For each value of P simulated, each Table contains: 
(i) the best fit parameters a and 6; (ii) the chi-squared for the fit X 2 %v an d (hi) the probability 
P(.X 2 > xju)- A s noted above, the closer the latter probability is to 1, the more consistent 
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the data-set is with the fitting function. A discussion of these results is given in the following 
Section. 



4 Discussion 

Ref. [4] examined the consequences of modifying the original linearly interpolating QuAdS 
Hamiltonian H(t) in eq. (00) to 

H(t) = H{t) + 6H(t) , (5) 

where the new term STC(t) has the form 

5H(t) = e(t)H E , (6) 

and the envelope function e(t) is required to vanish at t — and T. Ref. [I] gave 3 proposals 
for He, though it will not be necessary to review them here as we have a specific form in 
mind for STC(t) (see below). Eq. (JHJ) specifies a path in the space of 2 N x 2 N Hermitian 
matrices that begins and ends at Hi and Hp, respectively, and which by construction, differs 
from the linearly interpolating path specified by H(t). As noted in Section [TJ Ref. [I] showed 
that by doing such a path variation, a failure of QuAdS could be converted into a success. 

Including a noise interaction in the QuAdS Hamiltonian also causes the Hamiltonian 
path to deviate from the linearly interpolating path traced out by H(t). For non-uniform 
y-polarized noise, the noise term in eq. (j2J) is 

N 

8n{t) = -Y,mx , (7) 

where a y is the Pauli matrix o y for qubit j and Nj (t) y is the noise field that interacts with 
this qubit. Unlike the envelope function e(t), the noise fields {Nj(t) : j = 1, . . . ,N} need 
not vanish at t = or T. We now show, however, that for the noise used in the simulations 
presented here and in Ref. [8], the probability that a fluctuation is present at these times is 
small. To see this, note that for a fluctuation i to be present at t — (T), the fluctuation 
center ti must occur within a time r of t — (T) since the temporal width of the fluctuation 
is 2t. For our noise model, the time £j has a uniform probability distribution over the time 
interval [0, T] so the probability that ti is within r of t — (T) is r/T. Since n is the average 
fluctuation rate, the average number of fluctuations N T that occur in a time interval r is 
N T = nr. Since n = P/2a 2 r (see Section [2]), the total probability that a fluctuation is 
present at t = (T) is 

For P = 0.001 (0.013), the midpoint for the range of (T(N)) values found in the simulation 
is approximately 20 (39). Using this value for T in eq. (|SJ), and recalling that our simulations 
used a = 0.2 and r = 1 gives P tot = 6.25 x 10~ 4 (4.17 x 10~ 3 ) for P = 0.001 (0.013). Thus, 
with high probability, our noise interaction vanishes at t = (T), and our TC(t) is equal to 
Hi (Hp) at this time. 
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For non-uniform ^/-polarized noise, S7i(t) has N noise fields Nj(t), where 7 < A^ < 16 in 
our simulations. By comparison, the simulations in Ref. [8] used uniform noise in which all 
qubits see the same noise field Nj(t) = N(t), with j = 1, . . . , N. Thus our noise interaction 
with non-uniform y-polarized noise has an order of magnitude more variation parameters 
than the uniform noise used in Ref. [8j. Furthermore, since the non-uniform noise simulations 
were done to larger values of P than the simulations with uniform noise, each noise field in 
the former case has a larger average number of fluctuations Nf than in the latter case 
since Nf = PT/2a 2 r. Thus, because the simulations in the present paper were done using: 
(i) non-uniform noise; (ii) larger average noise power; and (iii) larger number of qubits, they 
contain larger Hamiltonian path variations than was possible in Ref. [8]. 

The simulation results presented in Tables [1] and [2] show that the scaling exponent b 
begins to decrease for P > 0.005. (As will be discussed below, we anticipate that there will 
be a maximum P value beyond which noise will begin to compromise QuAdS perfromance.) 
A second observation is that power-law scaling provides an excellent fit for all values of P 
simulated, while the exponential fit is not quite as good for P > 0.007. To compare our 
results with those in Ref. [8J a restricted power-law fit to the non-uniform noise results was 
done for 7 < JV < 12 and P = (0.001,0.003,0.005). This corresponds to the range of N 
and P values simulated in Ref. [8J for uniform ^/-polarized noise. The parameters for the 
restricted fit, together with the corresponding fit parameters for uniform ^/-polarized noise 
appear in Table [3J Note that a similar comparison is possible with exponential fits, though 
nothing new is learned and so we do not include that comparison here. From Table [3] we 
see that the scaling exponent is comparable for the two types of noise, with smaller b- values 
for non-uniform noise. Compared to the P = 0.000 results, we see that noise slows down 
QuAdS, but Table [3] gives a first indication that non-uniform noise may allow conditions for 
the slowdown to be ameliorated. 

A look at Figures [3] and |H shows that the initial rate of growth of (T(N)) appears to 
flatten out at intermediate values of N. This flattening out is less pronounced for P = 0.005, 
occurring over the range 12 < < 14; and is broader for P = 0.013, occurring over the 
range 10 — 11 < AT < 14. To test this observation we did separate fits for each of the data-sets 
in these Figures at small N (7 < N < 10) and (relatively) large N (13 < N < 16). The 
parameters for the two fits appear in Table HI As above, we only show results for a power-law 
fit. Table H] indicates that the small N fit grows at a faster rate (viz. larger b) than the large 
N fit. We see that the flattening out of (T(N)) at intermediate N marks a crossover from 
rapid initial growth to a region of slower growth. In an effort to further highlight this point, 
Figure EJreplots the data for P = 0.009, including the fits for small and large N. One clearly 
sees the data initially following the faster rising fit (b = 3.560) and then crossing over to the 
slower rising fit (b = 2.774). Although noise is plainly causing QuAdS to slowdown relative 
to noiseless QuAdS (see Figured]), Table H] indicates that non- uniform y-polarized noise can, 
for appropriate values of P and N, ameliorate the slowdown. One might wonder whether 
there is a connection between this ameliorating effect and the suggestion made in Ref. [4] that 
varying the Hamiltonian path away from the linearly interpolating path used in Refs. [2] and 
[3] might improve QuAdS performance. From that perspective, one might wonder whether, 
for 12 < N < 16, non-uniform y-polarized noise with 0.005 < P < 0.013 is inducing 
sufficient variation of the linearly interpolating Hamiltonian path to yield alternative paths 
with slightly larger minimum energy gaps, causing a (slightly) improved adiabaticity, and 
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so reducing (slightly) the noise-induced slowdown of QuAdS seen at smaller N values. A 
proper examination of this ansatz requires that the minimum gaps be determined for the USA 
instances and noise realizations that we simulated in this paper. We plan to carry out this 
analysis in future work. Note that for sufficiently large average noise power, noise-induced 
decoherence should ultimately rob QuAdS of its quantum performance-enhancements since it 
will cause the dynamics to crossover from quantum to classical. A quantitative determination 
of how this loss of quantum performance occurs presents a significant (though fascinating) 
challenge for future simulations. 
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Figure Captions 

Figure 1: Noiseless QuAdS simulation results for the median runtime T(N) (dimensionless 
units) versus number of bits N. The solid line is the best power-law fit to the data 
and the dash-dot line is the best exponential fit. The error bars give 95% confidence 
limits on the median. 

Figure 2: Simulation results for the noise-averaged median runtime (T(N)) (dimensionless 
units) versus number of bits N. The noise is polarized along y and has average noise 
power 0.001 and 0.003 (dimensionless units). Each datapoint is the median of 750 
runtimes (75 USA instances and 10 noise environments per USA instance). The solid- 
line is the best power-law fit to the data and the dash-dot line is the best exponential 
fit. The error bars give 95% confidence limits for each median. 

Figure 3: Simulation results for the noise-averaged median runtime (T(N)) (dimensionless 
units) versus number of bits N. The noise is polarized along y and has average noise 
power 0.005 and 0.007 (dimensionless units). Each datapoint is the median of 750 
runtimes (75 USA instances and 10 noise environments per USA instance). The solid- 
line is the best power-law fit to the data and the dash-dot line is the best exponential 
fit. The error bars give 95% confidence limits for each median. 

Figure 4: Simluation results for the noise-averaged median runtime (T(N)) (dimensionless 
units) versus number of bits N. The noise is polarized along y and has average noise 
power 0.009 and 0.013 (dimensionless units). Each datapoint is the median of 750 
runtimes (75 USA instances and 10 noise environments per USA instance). The solid- 
line is the best power-law fit to the data and the dash-dot line is the best exponential 
fit. The error bars give 95% confidence limits for each median. 

Figure 5: Re-plot of the simulation results for P = 0.009. The dashed (solid) curve is a 
restricted power-law fit (T(N)} = aN b through the data 7 < N < 10 (13 < N < 16). 
Similar plots are possible for the other P-values in Table HI though we do not include 
them to avoid repetition. 
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Table 1: Best fit parameters for power-law scaling (T(N)) = aN b for all values of average 
noise power P simulated. For comparison, best-fit parameters for noiseless QuAdS are also 
included. 



p 


a 


b 


x% 


P(x 2 > X%) 


0.000 


0.1016 


2.079 


0.321 


0.9999 


0.001 


0.07863 


2.200 


0.539 


0.9993 


0.003 


0.02595 


2.753 


0.869 


0.9967 


0.005 


0.02611 


2.816 


1.097 


0.9976 


0.007 


0.03783 


2.698 


2.919 


0.9393 


0.009 


0.07251 


2.458 


3.449 


0.9032 


0.013 


0.1636 


2.148 


2.689 


0.9523 



Table 2: Best fit parameters for exponential scaling (T(N)) = a[exp(bN) — 1] for all values of 
average noise power P simulated. For comparison, best-fit parameters for noiseless QuAdS 
are also included. 



p 


a 


b 


Xfu 


P{x 2 > X%) 


0.000 


4.707 


0.1282 


0.296 


0.9999 


0.001 


2.208 


0.1845 


0.185 


0.9999 


0.003 


1.580 


0.2302 


0.517 


0.9994 


0.005 


3.123 


0.1915 


3.307 


0.9136 


0.007 


3.924 


0.1799 


4.833 


0.7753 


0.009 


5.445 


0.1601 


4.927 


0.7653 


0.013 


7.112 


0.1440 


3.172 


0.9231 



Table 3: Comparison of power-law fit (T(iV)) = aN b for non-uniform and uniform (Ref. [8]) 
y-polarized noise for 7 < iV < 12 and average noise power P = 0.001, 0.003, and 0.005. 



p 


non-uniform 


uniform 




a b 


a b 


0.001 
0.003 
0.005 


1.107 x 10" 1 2.044 
5.095 x 10~ 2 2.439 
1.677 x 10" 2 3.020 


9.677 x 10" 2 2.108 
2.982 x 10~ 2 2.679 
6.607 x 10~ 3 3.429 
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Table 4: Restricted power-law fits (T(N)) = aN b for the P = 0.005,0.007,0.009,0.013 
data-sets for: (i) 7 < iV < 10; and (ii) 13 < iV < 16. 



p 


7 < N < 10 


13 < iV < 16 




a b 


a b 


0.005 
0.007 
0.009 
0.013 


3.064 x 10~ 2 2.730 
7.637 x 10~ 3 3.450 
7.086 x 10- 3 3.560 
3.201 x 10~ 2 2.924 


6.554 x 10" 2 2.465 
2.984 x lO" 2 2.767 
2.934 x 10~ 2 2.774 
3.060 x 10~ 2 2.759 
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Figure 1: 
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(a) P = 0.001 




(b) P = 0.003 



Figure 2: 
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(a) P = 0.009 




(b) P = 0.013 



Figure 4: 
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